STIMULATED RAMAN ADIABATIC PASSAGE 
IN THE FIELD OF FINITE DURATION PULSES 
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The theory of stimulated Raman adiabatic passage in a three-level A-scheme of the interaction of 
an atom or molecule with light, which takes the nonadiabatic processes at the beginning and the 
end of light pulses into account, is developed. 
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I. INTRODUCTION 



Adiabatic processes in atomic physics, owing to their 
stability with respect to a variation of parameters that 
describe the interaction with a field, play an especial role 
as a tool for manipulating atomic and molecular states. 
For instance, a fast adiabatic passage of the light pulse 
carrier frequency through the resonance with the atomic 
transition frequency allows one to obtain an atom or 
molecule in the excited state with a probability close to 
1 [H In this work, we study the stimulated Raman 
adiabatic passage (STIRAP), which can be realized in a 
three- and a multilevel scheme of the interaction between 
the atom and the field and was predicted as early as in 
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Figure 1: Scheme of the atom-field interaction. A Stokes 
pulse with the carrier frequency u>s affects the atom firstly 
and is followed by the second pumping pulse with the carrier 
frequency up, by partially overlapping in time with the first 
pulse. Here, W\ , Wi , and Wz are the energies of atomic states 
|1), |2), and |3), respectively 
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the 1980s d, El- The physical aspects and numerous ap- 
plications of STIRAP in various branches of physics and 
chemistry are discussed in reviews 0] and in work 0] • 
Recently, STIRAP has been demonstrated to occur in 
solids 0| ■ 

STIRAP is based on the existence of a trapped or 
"dark" state which arises provided that there is a two- 
photon resonance between an atom (in what follows, 
when speaking about an atom, we also mean a molecule) 
and a radiation field produced by two lasers [1, [t| [Tol ] . 
In the simplest case considered in this work, the laser 
carrier frequencies are so selected that they are close to 
the frequency of the transition between the excited state 
and two states with lower energies-both metastable or 
stable and metastable ones (the three-level A-scheme of 
the interaction between the atom and laser radiation). 
The radiation of one of the lasers - the pumping field - 
couples states |1) and |2) (see Fig. 1). The field of the 
other laser - the Stokes field - couples states |3) and |2). 
The difference between the frequencies of those fields co- 
incides with the transition frequency |1) <-> |3). In this 
case, some time after the interaction of the atom with the 
field has started, the probability to find the atom in the 
excited state |2) is close to zero, i.e. the atomic state is 
described by a linear superposition of the basic, |1), and 
metastable, |3), states. In this case, the populations of 
the states |1) and |3) are determined by the ratio between 
the intensities of Stokes and pumping fields. Being sub- 
jected to the action of the pumping field only, the atom 
is in state |3), whereas if only the Stokes field acts upon 
the atom, it is in state |1). If the intensity ratio changes 
slowly, the atom transits either from state |1) into state 
1 3) or vice versa. 

We adopt that, before the interaction between the 
atom and the field has started, the former is in state 
In this case, in order to transfer the population from 
state |1) to state |3) , it is necessary that a "counterintu- 
itive" sequence of pulses should affect the atom. Namely, 
the atom is first affected by the Stokes pulse and then by 
the pumping one which partially overlaps in time with 
the Stokes pulse. It is essential that, during the whole 
time interval with the atom-field interaction, the popula- 
tion of state 1 2) is very small, and the population losses 
owing to the spontaneous radiation emission from this 
state are insignificant. As a result, STIRAP provides the 
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population transfer between selected levels with a prob- 
ability close to 1. In addition, owing to the adiabaticity 
of the process (a slow variation of atom-field interaction 
parameters), the probability of population transfer prac- 
tically does not depend on wide-range variations of the 
shape and the intensity of light pulses. 

It is natural that the study of the influence of a nona- 
diabaticity, which reduces the population transfer prob- 
ability, on STIRAP has always drawn attention of re- 
searchers. First of all, it should be noted that the adia- 
baticity criterion for the process of interaction between 
the atom and the field was analyzed in practically every 
work, where STIRAP was considered (see, e.g., works 
0, II S, @1)- The criterion is based on the requirement 
that the eigenvectors of a Hamiltonian that describes the 
atom-field interaction should vary slowly in comparison 
with the difference between its eigenvalues Numer- 
ically, the adiabaticity is characterized by the parameter 
e which is reciprocal to the light pulse area. With the 
reduction of the adiabaticity parameter e, the population 
71.3(7;) of the target state |3) grows, in general, at t — > oo 
in the course of STIRAP (small oscillations are possible 
with a frequency of the order of the Rabi frequency of 
light pulses). 

In the most complete way, the dependence of 713(00) 
on e can be monitored in some cases where the shape of a 
pulse envelope allows an exact expression for the popula- 
tion of atomic states to be found [HI, [l3| . Provided that 
there are no losses in the atomic state population through 
the spontaneous radiation emission - this condition was 
postulated in the works cited above - the difference of 
the population in the target state |3) from 1 tends to 
zero, with a reduction of e, following different laws for 
different shapes of light pulses. For instance, if the am- 
plitudes of the Stokes pulse at t — > —00 and the pumping 
one at t — > 00 do not vanish, and those pulses can be de- 
scribed by analytical functions within the time interval 
[—00, 00] [HI, then 1 — 713 ~ exp (— /i/e), where /i is 
a certain constant of the order of 1. The dependence of 
such a type for the given class of functions has a general 
character, and the theory developed by Dykhne [l5|] and 
Davis and Pechukas [l6[ , which describes the law, follow- 
ing which the system tends to the adiabatic state with 
the reduction of e, is applicable to those functions. An- 
other type of dependence - a power-law with 1 — 77,3 ~ e 2 
- was found for pulses with a special shape, which al- 
lows an exact solution of the Schrodinger equation with 
nonzero first derivatives of the pumping pulse field at its 
beginning and of the Stokes pulse field at the moment 
of its termination to be obtained [H, HH ■ This result is 
based on a discontinuity of the derivative of the field at 
the time moment of its switching-on 

.13] 

, that is charac- 
teristic of other quantum-mechanical systems as well, in 
which the almost adiabatic evolution of the wave func- 
tion is possible [l7], EM HH ■ In the general case of the 
nonzero 77-th derivative at the beginning of a light pulse, 
taking the results of the cited works into account, one 
may expect that 1 — ti 3 ~ e 2n at e — > in the STIRAP 



case. 

While considering finite-duration pulses, we proceed 
from the fact that it is the only class of pulses which 
can be realized under real experimental conditions. We 
will examine the cases where the pulse damping can be 
neglected (short light pulses) and when the light pulse 
duration r considerably exceeds the inverse lifetime of 
an atom in the excited state, 7 . The former case was 
already analyzed for light pulses with identical ampli- 
tudes and with a special pulse shape that allowed the 
Schrodinger equation to be solved analytically [HI, [l3l ]. 
The latter case was analyzed earlier without making al- 
lowance for a nonadiabaticity associated with a jump of 
the field derivative at the beginning of the Stokes pulse 

MM- 

For short light pulses (7T <C 1) and the field-atom in- 
teraction close to the adiabatic one (e <C 1), we will find 
the target state population to an accuracy of e 2 for pulses 
with arbitrary shape, whose intensity grows in time at 
their beginning proportionally to t 2 , and to an accuracy 
of e 4 , if the intensity at the pulse beginning grows as 
t 4 . For long light pulses (jt 3> 1), we will demonstrate 
that the transient processes, which arise at the begin- 
ning of Stokes pulse action on the atom, though do not 
change considerably the probability of population trans- 
fer during STIRAP in comparison with the results ob- 
tained in works [2(J HH , do insert substantial corrections 
into them. 



II. BASIC EQUATIONS 

Consider an atom, the interaction of which with the 
field of two light pulses is described by the three-level 
scheme (Fig. 1). The pumping pulse denoted by the sub- 
script P below partially overlaps in time with the Stokes 
pulse (the subscript S). The Stokes pulse acts upon the 
atom firstly. The carrier frequency of the pumping pulse 
ujp is identical to that of the transition between states 
|1) and |2), and the carrier frequency of the Stokes pulse 
ujs to that of the transition between states |3) and |2): 

E = \'Ep{t)e- lUpt - ltpp{ ' t) + iE s (i)e-' ust - lw(l) +c.c. 

We suppose that the amplitudes of the Stokes, Ep(t), 
and pumping, E5 (t) , pulse fields change smoothly in time 
with a characteristic scale comparable to the pulse dura- 
tion T. 

State 1 1) , in which the atom stays before its interaction 
with the field, is considered to be stable or metastable, 
whereas state |3), into which we intend to transfer the 
atom using the STIRAP process, metastable, so that 
the variations of populations in those states within the 
time intervals comparable with the pulse duration, which 
take place owing to the processes of spontaneous emission 
from them, are neglected. Concerning the spontaneous 
emission from the excited state |2), we suppose that, in 
the course of this process, the atom transits into other 
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states different from |1) and |3) at the rate 7 (the life- 
time of an atom in the excited state is r sp = 7 -1 ). In 
this case, the atomic state can be described by a wave 
function, the time evolution of which is described by the 
Schrodinger equation 



if 1 



dt 



H(t)*(t) 



(1) 



with a non-Hermitian Hamiltonian that contains 7 
For resonant interaction between the atom and the field, 

uj p = (W 2 -Wi)/h, u s = (W 2 -W 3 )/H. 

The matrix representation of the atom-field interaction 
Hamiltonian in the rotating-wave approximation [l], Q 
and in the dipole approximation for the electric field 
strength looks like 



H(t) 



fl P (t) 

n P (t) -»7 n s (t) 
n s (t) 



(2) 



and the state vector is a column of probability am- 
plitudes c n (t) to find the atom in the state ipn{t) — 
exp(-iW n t/h) \n) (n = 1,2,3): 



V = [c 1 (t),c 2 {t),c 3 {t)} 1 



(3) 



Here, (l P (t) = -d l2 E P (t)/h and Q s (t) = -d 32 F, s (t)/n 
are the Rabi frequencies of the pumping and Stokes 
pulses, respectively; and d is the operator of atomic 
dipole moment. Without any loss of generality, the Rabi 
frequencies are considered to be real- valued [l(. We also 
assume that, at the moment when the pumping pulse 
starts to affect the atom, the latter is in state |1), i.e. the 
initial conditions look like 



ci{U 



1, 



c 2 (U 



0. 



c 3 (U 



0. 



(4) 



It is convenient to characterize the variation of the ra- 
tio between the Rabi frequencies of the Stokes and pump- 
ing pulses in time by the time dependence of the mixing 
angle d i] 



i)(t) = arctan 



n P {t) 



n s (t)' 

Then, Hamiltonian © looks like 



(5) 



H(0 = f 



n(i)cosi? 

fi(i)costf(i) -ry n(t)sini?(*) 
fi(i)sin#(i) 



where 



(0) 



(7) 



The adiabaticity parameter e, the reduction of which 
corresponds to the approach of the atom-field interac- 
tion to the adiabatic one, can be estimated as e ~ 



III. TIME DEPENDENCES OF LIGHT PULSE 
ENVELOPES 

We illustrate the accuracy of the results obtained be- 
low, which describe the dependence of the population of 
state 1 3) on the parameters of light pulses by compar- 
ing them with the results of numerical calculations for 
pulses with model shapes. Let the time dependence of 
the pumping pulse repeat that of the Stokes one, but 
with the delay td- 



n P (t) = n PQ F n (t-t d /2), 

n s (t)=Q S0 F n (t + t d /2), 



(8) 



where n = 1,2,3. 
tions 



enumerates the sequence of the func- 



F n (t) 



cos"(7rt/r), 
0, 



if \t\ < r/2; 
if |*| > r/2. 



(9) 



For such light pulses, the atom interacts with the field 
created by both pulses from the time moment U = — (t — 
td)/2 till the time moment tf = (r — td)/2, i.e. during the 
time interval t — td- If td > t, the Stokes and pumping 
pulses do not overlap in time. 

Among the whole family of light pulses ([8|), we use 
two envelopes, with n = 1 and n = 2. In the first case, 
the electric field strength has jumps of its first deriva- 
tive at the beginning and the end of pulses; in the lat- 
ter case, these are jumps of the second derivative. In 
addition, the first case for pulses with identical ampli- 
tudes, Slpo = flso = and a delay td = r/2 be- 
tween them is remarkable in that the Rabi frequency 
Q(t) — f2o does not depend on time during the time in- 
terval — t/4 < t < t/4, when both pulses interact with 
the atom simultaneously, and the mixing angle is a linear 
function of time, #(t) = irt/r + ir/4 [13J. This feature of 
the model makes it possible to obtain analytical expres- 
sions for integrals which appear in the theory and to use 
a simple example to illustrate the results obtained. The 
second model — it was applied earlier, e.g., in work [22j|— 
is close to Gaussian-like pulses which are often used for 
the simulation of a light pulse shape in theoretical calcu- 
lations @, SEE S3- 



IV. STIMULATED RAMAN ADIABATIC 
PASSAGE IN THE FIELD OF SHORT LIGHT 
PULSES 

In the case of short light pulses, the duration of which 
satisfies the condition jt <C 1, the term in Hamiltonian 
© which describes the relaxation can be neglected. Let 
us pass to the basis of characteristic (adiabatic) states of 
Hamiltonian ([6]). They satisfy the equation 



(10) 
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Simple calculations bring us to the eigenstates 
/2 

b (1 '- ] = ^-Oisintf- ^2 + V^cosi?) , (11) 



&( 1,0 ) = ip! cos i? — ip3 sin i?, 



(12) 



the matrix representation, the state vector looks like 
[ai._ (t), ai t o(t), ai j +(i)] T , and the Hamiltonian like 



H< l >(t) = £ 



-O iv^i? 
i%/2?9 ft 



(16) 



/2 

= — (-01 sin + ^ + V>3 cos t?) 
and the corresponding eigenvalues of Hamiltonian 



ti\i t ± — ± — hQ, 



h\ lfi = 0. 



(13) 



(14) 



Hereafter, to make notations short, we do not indicate 
the dependences of i9, ft, b^\ Aij (j — 0, ±), and the 
vectors of the rotating basis ip n (n = 1, 2, 3) on time (of 
course, if it does not cause misunderstanding). Index 1 
in the notations 6 < - lj - ) and Aij is introduced for the con- 
venience of subsequent calculations. It means the order 
of the adiabatic basis [further, we consider the adiabatic 
bases of higher (second and third) orders] . 

The state M 1 ' ) is known from the literature as a "dark" 
one, because an atom does not emit light from it [tl [20j. 

Provided that the sequence of pulses is "counterintu- 
itive" , i.e. when the atom is first subjected to the action 
of the Stokes pulse and then, with a certain delay, of 
the pumping one which continues to affect the atom for 
some time after the Stokes pulse terminates, the angle $ 
changes, according to formula (0, from zero to w/2. Tak- 
ing the initial conditions into account, we see that the 
atom is in the adiabatic state M 1,0 ) at the beginning of 
its interaction with the pumping field. If the parameters 
of light pulses change slowly enough, the atom remains 
in this adiabatic state during the whole period of the 
simultaneous interaction with the fields of both pulses. 
At the moment tf, when the Stokes field is switched off, 
■d = tt/2. As is seen from expression (fT2"]) . frt 1 - ) coincides 
with the state ^3 in this case, i.e. the atom transits from 
state |1) into state |3) due to its interaction with the field. 
The probability of population transfer |1) — > |3) is close 
to 1, provided that the process of interaction between the 
atom and the field is close to the adiabatic one. The stay 
of an atom or a molecule in the state b^- 1 ' ' during the 
whole period of its interaction with the field composes 
the physical basis of STIRAP. 

Using the basis of adiabatic states, the wave function 
can be written down in the form 



(15) 



j=o,± 



Here, aij(t) is the probability amplitude of finding the 
atom in the j-th adiabatic state. Substituting func- 
tion (fT5|) into the Schrodinger equation (JXJ) , we obtain 
the Schrodinger equation in the adiabatic basis. In 



The nonadiabaticity of the atom-field interaction is de- 
scribed by non-diagonal elements of Hamiltonian Q16p . 
Taking into account that ~ 1/r, we see that non- 
diagonal elements are about eft by the order of mag- 
nitude. If they are neglected, the vector of atomic state 
in the adiabatic approximation looks as 



1 / Au(f)df • (17) 



*= E ai,j(*i)& (1 ' j) (*)exp( 

j=0,± 



One can see that, to within the phase, the amplitudes 
of adiabatic states remain constant during the whole pe- 
riod of the atom-field interaction. If the atom was in 
the "dark" state i/ 1,0 ) at the beginning of its interaction 
with the pumping pulse (the time moment ti), it stays 
in it after the interaction terminates. Hence, in the adi- 
abatic approximation, we have the population transfer 
between states |1) and |3) with the probability equal to 
1. 

Now, let us find small nonadiabaticity-induced correc- 
tions to the population transfer from state |1) to state |3) 
in the general case of an arbitrary pulse shape. For this 
purpose, it is necessary to take into consideration that 
the derivative of 1? in Hamiltonian (|16p differs from zero. 
Let us take advantage of the formalism used for the de- 
scription of the quantum-mechanical system in adiabatic 
bases of higher orders, which is well-known from the lit- 
erature (see [14], [2J, [25| )• Let us pass to the basis of 
eigenstates of Hamiltonian (fTo]) 



2ft 



ft 



ft-ft 

2ft 



b^+\ (18) 



6(2-0) = + «6d,o) + 4V^>, (19) 

ft n ft 



6 ( 2 ,+) = » - » b (i,-) , iVM b d,o) + » b a,+) (20) 
2fi ft 2C1 

with corresponding eigenvalues 

h\ 2t± = ±Ui, h\ 2 , = 0. (21) 

In Eqs. (|18[) - (|2Tj) , we introduced the notation 



(22) 



ft = \/ft 2 + Ad 2 . 
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In the basis of adiabatic states the wave function 

can be written down in the form analogous to expression 
([TBI with the substitution 1 — > 2 



i=o,± 



(23) 



Here, a 2 j(t) is the probability amplitude of finding the 
atom in the state b^\ In the matrix representation, 
the state vector in the basis of states b^ 2 '^ looks like 
[a 2> - (i), O2,o(*), 02,+ (i)] T , and the Hamiltonian like 



H< 2 )(i) 



where the notation 
2V2 





-n 


i/3 





h 


-if] 





-i/3 


2 





if3 





(24) 



(25) 



was used. Neglecting the non-diagonal elements in 
Hamiltonian (f24|) - or, equivalently, the dependence of 
&2,j(t) on time, - we obtain the wave function of an atom 
in the form of a superposition of characteristic states of 
Hamiltonian (Til 



*= a2,j(U)b (2 ' j \t)exp(-i f \ 2d (t')dA. (26) 
i=o,± V I J 




40 60 80 100 120 140 160 180 200 

flpOT 

Figure 2: Dependences of the probability of population trans- 
fer n3 from the atomic state |1) into state |3) on the Rabi fre- 
quency of a pumping pulse Qpo measured in 1/r-units in the 
field of light pulses of form JS), (|9]) with n = 1 and td = r/2 
for various ratios between Qpo and Qso, calculated by formula 
(|33p and by the numerical integration of the Schrodinger equa- 
tion (O with Hamiltonian ([2]). For curve 1 with Qpo = Qso, 
both calculation methods give identical results. Curves 2 
were calculated for Qso = 2Qpo, curves 3 for Qso = 5Qpo- 
Dashed curves denote the results of numerical integration of 
the Schrodinger equation 



[y/2 

-C2,+ (ti) I — COST? 



jV2tfsintf 



i*(t). 



Passing from the basis b^ 2 ^\t) to the basis 6^^(i) and, 
then, to ^> n with the use of relations (fTTj) - (fT3"|) and (fT5|) - 
(I20p . we find the population amplitudes c n (i) for states 



ci(t) = o 2 ,_(*;) I — smt?+ ~ I e 49W - 



. . / y/2 . zV2tfcostf\ _ i4(t) 
-a 2 ,+ (*i) I —sin ~ je 



c 2 (t) = -aa.-CiO^e**^ + a 2)0 (^+ 



-Oa.+ fe^sinife-'*^ 



(27) 



(28) 



C3(*) = aa,_(*i) I -y cost? - ^— ^ 1 e 1 *^- 



-O2,o(*i) 



if2 sin 



where 



*(*) = 2 / fi W 



(29) 



(30) 



The cocfhcicnts a 2 j(ti) are determined from the initial 
conditions (fj| 



02,- = 02,+ (*i) = - 



fi(ti) ^ 



O2,o(*j 



fi(*<) ' 



(31) 
(32) 



After the Stokes pulse terminates, the population of 
state 1 3) does not change any more in time. From formu- 
las (J2ZH32J), we find that 



n 3 



dt) 



(33) 
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According to the result obtained, the probability 723 = 
|c 3 | 2 of population transfer |1) — * |3) is governed by the 
first derivative of the mixing angle ■& (formula ([5|)) at the 
beginning of a pumping pulse (time U) 



n s (t ( )' 

and at the end of a Stokes pulse (time tf) 



*(*/) 



(34) 



(35) 



From formulas 
maximal value, if 



fip(ti) =±(l S (tf). 



SJ), it follows that U3 reaches its 



(36) 



Condition ([36]) shows that the optimum conditions for 
population transfer are obtained, in particular, for sym- 
metric, with respect to the maximum, pulses with identi- 
cal shapes and amplitudes. In this case, provided that the 
pulse amplitude or the pulse delay is selected properly, 
so that the cosine in Eq. (j3"3"|) is equal to 1, the popula- 
tion transfer from state |1) into state |3) is complete in 
the approximation of adiabatic atomic evolution in the 
basis of states 6^ 2j '(t). At the same time, the deviation 
of 713 from 1 with the variation of the cosine argument 
can reach 16-d (U) 2 /ft(ti) 2 , which is of the order of e 2 . 

Note that, provided that the mixing angle is propor- 
tional to time during the simultaneous action of light 
pulses on the atom and the frequency f2 does not depend 
on time, the non-diagonal elements in Hamiltonian (|24p 
are equal to zero, and functions (f27|) -(l30 l) are the exact 
solutions of the Schrodinger equation [l3j. 

In Fig. 2, the results calculated by formula (|3"3"|) for 
pulses with envelope © with n = 1 and t& = r/2 are 
depicted for various ratios between fipo and flso- They 
are also compared with the results of numerical integra- 
tion of the Schrodinger equation |T]) with Hamiltonian @ 
and the same parameters of the atom-field interaction. 
For curve 1 corresponding to flpo — f^soi both calcula- 
tion methods give an identical result, because expression 
(I33p is the exact solution of the Schrodinger equation 
in this case. As is seen from the figure, the probability 
of population transfer tends to 1 with increase of Qpqt, 
which corresponds to a reduction of the adiabaticity pa- 
rameter e. Simultaneously, the discrepancy between the 
results of calculations by formula (|33[) and numerical in- 
tegration of the Schrodinger equation decreases. 

Equation (|33|) includes the first derivatives of the mix- 
ing angle ■& with respect to time calculated at the begin- 
ning of the pumping pulse and the end of the Stokes one. 
If the derivatives are equal to zero, then 713 = 1. In this 
case, in order to find an expression for 713, which would 
make allowance for the nonadiabaticity of the atom-field 
interaction, we should seek the wave function in the adi- 
abatic basis of higher order than that of b^ 2 ^\t). 

Now, let us analyze how the non-zero second deriva- 
tives of $ at the beginning of the pumping pulse and at 



the end of the Stokes one - either or both - affect the 
probability of population transfer between states |1) and 
1 3) under conditions close to those of the adiabatic atom- 
field interaction. The calculation routine is the same, as 
was used at the derivation of formula (|33[) for the popula- 
tion of state |3). Being interested only in the case where 
the proximity of the population transfer to the adiabatic 
one is governed by the second derivative of the pumping 
field at the beginning of the pumping pulse and the sec- 
ond derivative of the Stokes field at the end of the Stokes 
pulse, we assume that 



0(t<) = 0(t f ) = 0. 



(37) 



We obtain the eigenvalues fiX^j (j — ±, 0) and eigen- 
statcs M 3j )(i) of Hamiltonian and take the m as a 
new basis. In the Hamiltonian H^(t) which describes 
the atom-field interaction in this basis, we neglect non- 
diagonal elements. In this basis, the wave function looks 
like (|26[) with index 2 being substituted by 3. With 
regard for the relations between basis vectors b^'i'it), 
&( 2 -?')(i), b^' j \t), and ip n (n = 1, 2,3), as well as the ini- 
tial conditions ([3]), we find the population of state |3): 



n 3 



where 





n(u) 2 n(t f ) 2 + i6i?(ti)i?(t/)cos$ 


2 




(n(t,)< + le^/) 2 ) 



, (38) 



I = ~ j Jn 2 + m(m-my n~ 4 dt. (39) 



The second derivatives of the mixing angle at the be- 
ginning of the pumping pulse and at the end of the Stokes 
one are coupled with the second derivatives of the corre- 
sponding Rabi frequencies, 



n P (tj) 

Os(ii) : 



m) = - 



(40) 



From expression 
imal value, when 



it follows that n 3 reaches the max- 

n P {u) 2 3(t f ) = ±n s (t f ) 2 #(t i ). (4i) 



As is seen from formula (|41|) , the maximal population 
transfer between states |1) and |3) is achieved, in par- 
ticular, in the field of pulses, for which the values of 
are identical at the beginning of the pumping pulse and 
at the end of the Stokes one, = |i?(i/)|, and the 

Rabi frequencies are also identical at these time moments, 
f2(ti) = fi(t/). In this case, provided that the intensities 
of pulses are selected properly, so that $ = 27m, where n 
is an integer number, the probability of population trans- 
fer is close to 1 . The deviation of the probability from 1 
with the variation of $ is about e by an order of mag- 
nitude. 
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QpQT 

Figure 3: Dependences of the probability of population trans- 
fer n-i from atomic state [1) into state |3) on the Rabi fre- 
quency of a pumping pulse Qp measured in 1/r-units in the 
field of light pulses of form J8), © with n = 2 and td = r/2 
for various ratios between Qpo and Qso, calculated by for- 
mula (|38p and by the numerical integration of the Schrodinger 
equation (fTJ) with Hamiltonian For curves 1, Q.po = fiso; 
curves 2 were calculated for Q$o = 2fipo, an d curves 5 for 
fiso = 5f2po- Dashed curves denote the results of numerical 
integration of the Schrodinger equation 



transfer, and they give U3 = 1. To find these corrections, 
it is necessary to pass to the adiabatic basis of higher 
order. Should the first derivative different from zero at 
the moments U and tf be a derivative of the n-th order, 
the order of magnitude of the nonadiabaticity correction 
would be e 2n . 



V. STIMULATED RAMAN ADIABATIC 
PASSAGE IN THE FIELD OF LONG LIGHT 
PULSES 

Now, consider the population transfer between atomic 
states |1) and |3) in the course of STIRAP in the case of 
long light pulses, when the time of the atom-field inter- 
action considerably exceeds the time of the spontaneous 
emission from the excited state, 77 > 1. Consider the 
"bright", ipb, excited, ip e , and "dark", ipd, states denned 
by formulas [2(| 

ijj b = sin 0(t)V>i + cos 0(f)t/>3, (42) 



1p e =1p2, (43) 



In Fig. 3, the results of calculations by formula |38J for 
pulses of form ([S]),© with n = 2 and td — r/2 for vari- 
ous ratios between flpo and ilso are depicted and com- 
pared with the results obtained by numerical integration 
of the Schrodinger equation |T]) with Hamiltonian @ and 
the same parameters of the atom-field interaction. As is 
seen, the maximal values of 713 are reached for pulses 
with |i?(ii)| = and fi(£j) = fi(t/), in accordance 

with the analysis of expression (|38p given above for the 
probability of population transfer. The oscillations of the 
population expectedly decrease with increase of the pulse 
area, i.e. as the interaction with the field approaches the 
adiabatic one. 

The obtained results [formulas ((33)) and (|38j) ] describe 
the dependence of the population transfer on the light 
pulse parameters in the case where the interaction be- 
tween the atom and the field is close to the adiabatic one. 
As is seen from Figs. 2 and 3, they correctly describe the 
oscillating dependence of the population transfer proba- 
bility in this range of parameters. In Fig. 3 which corre- 
sponds to nonzero second derivatives of the mixing angle 
at the moments ti and t f , the plotted curves are located 
much closer to the value = 1, than the corresponding 
curves in Fig. 2, where already the first derivatives of the 
mixing angle are different from zero. It was to be ex- 
pected, because the amplitudes of oscillations are of the 
order of e 4 and e 2 for the dependences plotted in Figs. 3 
and 2, respectively. 

If both the first and second derivatives of the mix- 
ing angle equal zero at the time moments U and tf, ex- 
pressions (|33[) and (|38p do not contain any more nona- 
diabaticity corrections to the probability of population 



ipd = costf(*)V>i - sini?(t)V>3, 



(44) 



where ipj (j — 1,2,3) are the basis wave functions of 
the rotating basis, in which Hamiltonian J5J is written 
down. The functions ipj differ from the functions \ j) only 
by time-dependent phases. Only one of those states, ipd, 
which coincides with (jT2J) , is the eigenstate of Hamilto- 
nian ([2]). At the beginning of the atom-field interaction, 
$ = 0, and the atom is in the state ipd- Should it stay in 
this state during the whole period of the interaction with 
the field, then, at the moment, when the Stokes pulse 
terminates and $ grows up to z9 = 7r/2, the population 
would be completely transferred from the state ip\ into 
the state ip^. 

The vector of state constructed of the probability am- 
plitudes Cb, C e , and Cd to find the atom in the states 
ipb, ip e , and ipd, respectively, looks like [Cb, C e , Cd] T , and 
the Hamiltonian in this basis is 



H ( bcd )(t) 



n m 
n -«7 

-id 



(45) 



The probability of population transfer from the state ip± 
into the state 1P3 — or, equivalently, the population 71,3 of 
the state ips, which we are interested in - is equal to 



^3 



|C 3 (*/)r 



(46) 



If STIRAP occurs in the field of long pulses, the as- 
sumption is usually made that the characteristic time of 
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a probability amplitude variation is of the same order of 
magnitude as the duration of light pulses [13, Hi], |22| . 
This assumption allows the Schrodinger equation to 
be solved by the iteration method, supposing that the 
derivatives of amplitudes are small in comparison with 
ft. In essence, the transient processes that start at the 
beginning ti of the interaction between the atom and the 
pumping pulse are neglected. This approach always gives 
rise to the correct first term in the expansion of the prob- 
ability in a series in the adiabaticity parameter e. 

In order to take the correction to the population 
transfer probability associated with the transient pro- 
cesses arising at the switching-on of a pumping pulse 
into account, we solve the Schrodinger equation in two 
stages. First, we suppose that the left-hand side of the 
Schrodinger equation with Hamiltonian (I45f> represented 
in the basis of states ?pb, "0e, and ip d has the same or- 
der of magnitude as the term proportional to Q on the 
right-hand side, and solve the equation in the time in- 
terval [it,*i], where 7 _1 <C t\ — U <C r. Such a ii-value 
can always be found, bearing in mind the condition of 
long interaction between the atom and the field, jt ^> 1 . 
When solving the Schrodinger equation in this time inter- 
val, we neglect the time dependence Cl(t) and adopt that 
f2(t) = fi(tj). As a result, we take damped oscillations 
of the amplitudes at the beginning of the atom-field in- 
teraction into account. Further, we solve the Schrodinger 
equation in the interval \t-\ , t fl, by supposing now, as was 
done in works [13, HH, H2| , that the characteristic time of 
amplitude derivative variations has an order of magni- 
tude of T. 

Let us pass to the variables 

rjd = hiC d , Vb = C b /C d , r) b = C b /C d . (47) 

From the Schrodinger equation Jl} with Hamiltonian 
(gnj), we find 



f]b = --zSlrje - r] b f] d + 



(48) 



(at the end of calculations, we put e = 1) and seek rjb, 
r/ e , and r\ d in the form 



ri b = ^2Hb, n e n , 

71=0 



n=0 



Vd = ^2 H d . n e r 



(52) 



(53) 



(54) 



First, let us consider the case f}(ti) ^ 0, i.e. when 
the Rabi frequency of a pumping pulse at its beginning 
linearly depends on time (see Eq. (f3~4"|) ). Substituting 
Eqs. ([32]) -([54]) in Eqs. (|48]h([50]) with regard for the ini- 
tial conditions 



rib = rje = Vb = 0, 



(55) 



which follow from Eq. ([4]), we find, after simple calcula- 
tions, that 



H d .o — H d ,i — 0, 



(56) 



27a 2 2a 2 7 1 ,/ ft' rz — - \ 



* 2 (7 2 



20?) e-ht' 



a 2 Ja 2 



1^2 

4 ' 



•Mill 7p/ fi < 



h 2 



(57) 



which is necessary for further calculations of the popula- 
tion in the state ip d . Here, the notations 



n t = n(u), a = #{u), t' = t-u 



(58) 



*o 1 



fid = -Vbd- 



(49) 



(50) 



The population of state |3) after the Stokes pulse termi- 
nates coincides with the population in the state ip d , being 
equal to 



n 3 = exp [2r7d(i/)] = exp / 2f] d (t)dt. 



(51) 



are used. 

Now, consider the time interval [ti, £/], where oscilla- 
tions of the population in the states ipb, ipe, & n d ip d prac- 
tically disappear (since jti 3> 1), and the amplitudes 
of those states slowly change in time and approximately 
follow the variations of light pulse Rabi frequencies. Sim- 
ilarly to what was done when solving Eqs. (I48[) — (f5U|) in 
the time interval [ii,ti], we formally introduce a small 
parameter e into them to mark the magnitude of coef- 
ficients. Since max(f2r) 1 and 77 > 1, let us write 
down Eqs. I[i8] ) -([5D ]) in the form 



(59) 



fib = --Slr) e e 1 - fibf] d + 1?, 



Consider Eqs. (I48[) (f5U|) in the time interval [ij,ti], 
where f) is small. To mark this smallness, let us formally 
introduce the parameter e <C 1 at in those equations 



-Ctr) b e 1 



VeVd, 



(60) 
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Figure 4: Dependences of the probability of population trans- 
fer 713 from the atomic state |1) into state |3) on Q.qt in 
the field of light pulses of form p), (JH with n — 1 and 
t d = t/2 calculated for various jt = 10 (1), 20 (2), 40 (3), 
and 100 (4) by formula (|66[) and by the numerical integra- 
tion of the Schrodinger equation fl]) with Hamiltonian ([2j|. 
Thick solid curves are the results of numerical integration of 
the Schrodinger equation, thin solid curves are the results of 
calculations by formula (|66[1 . dashed curves are the results of 
calculations by the same formula without regard for the first 
term in the exponent which is responsible for a nonadiabatic- 
ity that is inserted by the jump of the first derivative of $ at 
the beginning of a pumping pulse 



Vd = -ribti- 



(61) 



The absence of the factor e _1 at r]b, Ve, and % means 
that the characteristic times of their variations are of the 
order of the light pulse length r. 

The solution of Eqs. ([59|h(l6l]) is sought in form (l52|) - 
(f54|) . After simple calculations, we obtain 



Hd.o — 0, 



Hd.i — 



2 7 tf 2 



H, 



d,2 



Ad 2 tl (n 2 - 2 7 2 ) Add ( 7 2 



Q 2 



n 5 



n 4 



(62) 



(63) 



The expression for i)d obtained from Eqs. (|54|) and ((56)) - 
([57)1 for ti<t<t\ has to transform at the time moment 
t = t\ into the corresponding expression obtained for 
ti <t<tf. Really, the quantity Hd,i(ii) from Eq. l|6"2")) 
is equal to Hd^{t{) from Eq. ((57)) . because the oscillating 
terms in the latter, due to the exponential damping, are 
practically zeroed at this moment. No misunderstand- 
ing should invoke the comparison made between terms 
of different e-orders, because this parameter is different 
at different time intervals: in the former case, it marks 
terms with the $ order of magnitude; in the latter, the 
terms of the order of the adiabaticity parameter e. 



The terms that correspond to H^p from Eq. (|63|) and 
have higher e-orders than those presented in Eqs. (|56")) - 



(|57f have to appear in the interval [ti, t\]. The indicated 
order is the maximal one in the expansion of fjd that 
should be taken into account for a linear time dependence 
of the Rabi frequency of the pumping field, because the 
exceeding of the calculation precision occurs otherwise. 

Now, let us calculate the probability of population 
transfer from state |1) into state |3). From Eqs. (I5ip . 
([57). ([52]). and ([M]), we find 



ri3 = exp 



8a 2 , , 



4 7 



d 2 



v 2 h (n 2 ~ 2 7 2 ) w ( 7 2 - n 2 ) 



O 4 



dt 



(64) 



Here, the first term in the exponent emerges owing to 
damped oscillations of the population with a frequency 
of the order of f2i which arise at the beginning of a pump- 
ing pulse. The other terms obtained in works [20L l2ll] are 
associated with the quasistationary evolution of the pop- 
ulations of atomic states with a characteristic time of the 
order of the light pulse duration. 

Let us illustrate the obtained result in the case where 
the integrals in Eq. ([64)) can be calculated analytically. 
Consider light pulses of form ([8]) , ([9]) with n = 1 and at 
flpQ = ttso = f^o and td = t/2. For such pulses, 



d 



-t. 



(65) 



and the second integral in Eq. ()64j) vanishes, because = 
and 1? = 0. Simple calculations bring about 



ri3 = exp 



8tt 2 



(7 2 



n 2 ) 



2 7 7T 2 



(66) 



As is seen, the relative contribution of expression (|66|) 
to the exponent, which is associated with transient pro- 
cesses at the beginning of a pumping pulse, is of the or- 
der of 4/( 7 r). For example, at 7 r = 40, in the case of 
the atom-field interaction close to the adiabatic one, i.e. 
QoT 1, the corresponding correction to the quantity 
1 — «3 is about 10%. 

In Fig. 4, the results of numerical calculations of the 
probability of population transfer from state |1) into 
state 1 3) obtained by the numerical integration of the 
Schrodinger equation are shown, as well as the results 
of calculations by formula (|66p. where allowance is 
made or not for the first term in the exponent which 
is responsible for a nonadiabaticity inserted by the 
jump in the first derivative of the mixing angle -d at the 
beginning of a pumping pulse. The figure demonstrates 
that taking the nonadiabaticity associated with the 
jump of at the time moment ti into consideration 
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substantially improves the accuracy of calculations. 
The dependence of the population transfer probability 
on the light pulse area obtained in such a way (in Fig. 4, 
the light pulse area is parametrized by the product 
Qot) practically coincides with the result of numerical 
calculations of this quantity from the Schrodinger equa- 
tion at 7 > 20/r. The result of calculation by formula 
(prJ)) at 7T = 10 reproduces - though on the average 

- the corresponding result obtained by the numerical 
integration of the Schrodinger equation. At the same 
time, they appreciably differ from each other, because 
the condition that there must be such t\, which satisfies 
both inequalities -^{ti — U) 3> 1 and {t\ — U) -C (£/ — *i) 

- simultaneously, is violated. Really, in the case of the 
light pulses under consideration, (tf — ti) — r/2, so that 
7(ii — ti) cannot exceed 5. As a result, as is seen from 
Eq. (|57|) , the exponent that is responsible for the damp- 
ing of population amplitude oscillations in the dark state 
amounts to only e -125 at the end of the simultaneous 
interaction of the atom with the fields of Stokes and 
pumping pulses, which contradicts the assumption about 
the oscillation termination within a short, in comparison 
with r/2, time of the atom-field interaction, which is 
necessary for expression (|66[) to be valid. 

In the case where the first derivative of the mixing an- 
gle i9 at the beginning ti of a pumping pulse is equal to 
zero, whereas the second derivative is nonzero, it is also 
possible, within the calculation scheme described above, 
to obtain an expression for the population transfer prob- 
ability similar to formula (|64[) . The correction in the 
exponent, which emerges due to damped population os- 
cillations arising at the beginning of a pumping pulse, is 
of the order of e 4 in this case [in expression (|64[) . it is 
of the order of e 2 ], which, in general, is much less than 
the values of integrals included into Eq. ([M]) . Hence, it 
is eligible to neglect the transient processes arising at the 
beginning of a pumping pulse, provided that the time 



dependence of the mixing angle is described by a power 
law i? ~ t n with n > 2. In this case, the probability 
of population transfer can be found in the quasistation- 
ary approximation, at least with an accuracy of not worse 
than e 3 , supposing that the characteristic variation times 
of atomic state populations are of the order of the light 
pulse duration [2(| Hl| . 

VI. CONCLUSIONS 

We have analyzed the influence of the extra nonadi- 
abaticity associated with the non-analytical behavior of 
the field strengths of light pulses at the beginning and 
the end of their action upon the atom in the course of 
stimulated Raman adiabatic passage on the probability 
of population transfer. The cases where the light pulses 
are much shorter than the lifetime of an atom in the 
intermediate state and when the time of the atom-field 
interaction considerably exceeds the time of the sponta- 
neous emission in this state, have been considered. In 
both cases, the additional nonadiabaticity is maximal, 
when the field strength grows linearly (or the intensity 
quadratically) with time at the beginning of a pumping 
pulse. For short light pulses, the optimum conditions for 
population transfer are reached, if the time-derivative of 
the Rabi frequency of a pumping pulse at the time mo- 
ment of switching-on is equal to that of a Stokes pulse 
at the time moment of switching-off. In the case of long 
light pulses, the correction, which is related to the tran- 
sient processes occurring at the beginning of a pumping 
pulse, to the theory developed in works [2(1 HlJ can ap- 
preciably change the result only if the first derivative of 
the pumping field strength differs from zero. 

The work was executed in the framework of the themes 
Nos. V136 and VTsl39 of the NAS of Ukraine. 
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